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ABSTRACT 

Recent studies have shown that strong correlations are observed between the low frequencies (1-10 Hz) of 
quasi-periodic oscillations (QPOs) and the spectral power law index of several black hole (BH) candidate sources, 
Qi in low (hard) states, steep power law (soft) states, and transitions between these states. The observations indicate 

that the X-ray spectra of such state (phases) show the presence of a power-law component and are sometimes 
related to simultaneous radio emission, indicating the probable presence of a jet. Strong QPOs (>20% rms) are 
present in the power density spectrum in the spectral range where the power-law component is dominant (i.e., 
60%-90%). This evidence contradicts the dominant, long-standing interpretation of QPOs as a signature of the 
thermal accretion disk. We present the data from the literature and our own data to illustrate the dominance of 
power-law index-QPO frequency correlations. We provide a model that identifies and explains the origin of the 
QPOs and how they are imprinted on the properties of the power-law flux component. We argue for the existence 
of a bounded compact coronal region that is a natural consequence of the adjustment of the Keplerian disk flow to 
the innermost sub-Keplerian boundary conditions near the central object and that ultimately leads to the formation 
of a transition layer (TL) between the adjustment radius and the innermost boundary. The model predicts two 
phases or states dictated by the photon upscattering produced in the TL: (1) a hard state, in which the TL is 
optically thin and very hot (kT £ 50 keV), producing photon upscattering via thermal Comptonization (the photon 
spectrum index T~1.7 for this state is dictated by gravitational energy release and Compton cooling in an 
optically thin shock near the adjustment radius), and (2) a soft state that is optically thick and relatively cold 
(kT < 5 keV; the index for this state, T ~ 2.8, is determined by soft-photon upscattering and photon trapping in a 
converging flow into the BH). In the TL model for the corona, the QPO frequency i/ high is related to the 
gravitational (close to Keplerian) frequency i/ K at the outer (adjustment) radius and i/ low is related to the TL’s 
normal mode (magnetoacoustic) oscillation frequency i^ma. The observed correlations between index and low and 
high QPO frequencies are readily explained in terms of this model. We also suggest a new method for evaluation 
of the BH mass using the index-frequency correlation. 

Subject headings: black hole physics — stars: oscillations — X-rays: general 


1. INTRODUCTION 

The short dynamical timescales and the coherence of quasi- 
periodic oscillations (QPOs) observed in accreting X-ray 
bursts strongly suggest that QPOs can provide valuable in- 
formation on the accretion dynamics in the innermost parts of 
these systems. In this respect the discovery of low 20-50 Hz 
QPOs in luminous neutron star (NS) binaries by van der Klis 
et al. (1985), kilohertz QPOs in NSs by Strohmayer et al. 
(1996), and hectohertz QPOs in black holes (BHs) by Morgan 
et al. (1997) opened anew era in the study of the dynamics near 
compact objects. Following these discoveries Psaltis et al. 
(1999) demonstrated that these NS and BH low and high fre- 
quencies follow a remarkably tight correlation. These features 
are the horizontal-branch oscillations (HBOs), along with the 
low-frequency noise Lorentzian j/ 1ow and the lower kHz QPO 
^high f° r NS and BH, respectively. 
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Belloni et al. (2002) have updated the correlation of Psaltis 
et al. (1999) adding data from Nowak (2000), Boirin et al. 
(2000), Homan et al. (2001), Di Salvo et al. (2001), and 
Nowak et al. (2002). Psaltis et al. suggest that the low and 
high frequencies correlate in a way that seems to depend only 
weakly on the properties of the sources, such as the mass, the 
magnetic field, or possibly the presence of a hard surface in 
compact objects. Mauche (2002) has reported low- and high- 
frequency QPOs in the dwarf nova SS Cygni and in VW Hyi 
(see Wouldt & Warner 2002 for the VW Hyi observations) and 
has called attention to the fact that these QPOs extend the 
correlation of Psaltis et al. downward in frequency by more 
than 2 orders of magnitude. This frequency correlation ob- 
served over such a broad range of sources indicates that a 
common phenomenon is responsible for QPOs in compact 
objects and raises the question of whether one can find an 
observational relation that would distinguish BH candidates 
from other compact sources. 

The comparative study of spectral and QPO features in BH 
candidate sources has yielded important information that can 
help answer this question. As far as spectral properties are 
concerned, there is increasing observational evidence that the 
large complex of spectral states originally devised by Belloni 
et al. (2000) to classify the highly variable X-ray spectra of BH 
Galactic sources like GRS 1915+105 can be reduced to a few 
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simple configurations (see, e.g., McClintock & Remillard 2003; 
Belloni 2003): (1) a low, hard power law-dominated state in 
which the photon spectral index is near 1.5; (2) a steep power 
law state in which the index is about 2.7; and (3) a thermally 
dominated state. The observed variability is explained by 
transitions between these three canonical configurations. 

Low-frequency (1-10 Hz) QPOs are generally observed to 
be present in states 1 and 2, during transitions between states 1 
and 2, but never in state 3. Both high- and low-frequency 
QPOs have been thought to be associated with the accretion 
disk, although their specific origin has not been widely un- 
derstood (but see Titarchuk & Wood 2002, hereafter TW02, 
for an explanation of their origin). Attempts to find consistent 
correlations with observable disk parameters have led to 
mixed results. In particular, there are so many exceptions to 
the observed correlations of frequency to disk spectral values 
sometimes observed that no general statement about the re- 
lation of the QPO with any spectral parameter solely associ- 
ated with the disk, such as thermal flux, temperature, inner 
radius, etc., can be made. For example, exceptions to one of 
the most widely quoted correlations, i.e., between disk flux 
and frequency, have been presented for the well-studied mi- 
croquasar GRS 1915+105 by Fiorito et al. (2003) and for XTE 
1550-564 by Remillard et al. (2002a, 2002b). 

On the other hand, a number of recent studies have revealed 
consistent and robust correlations between the photon index of 
the power-law component of the X-ray spectrum and low QPO 
frequency. Such strong correlations are observed the well- 
studied BH binary GRS 1915+105 in almost all situations in 
which the QPO is observable and in several other BHs over a 
wide range of observations and states. 

These studies include the following: 

1. A comprehensive study of correlations between photon 
index and QPO frequency in BH sources by Vignarca et al. 
(2003, hereafter V03); these results show strong correlations 
between index and QPO frequency, as well as saturation of the 
index between at the values 1.7 ±0.1 and 2.7 ± 0.2 for several 
BH sources and for the same source in different states. 

2. Recent studies by Kalemci (2002) for a number of Ga- 
lactic BH X-ray transients observed during outburst decay that 
confirm the persistence of the correlation between frequency 
and index during state transitions. 

3. Our own studies of GRS 1915+105 on x states (Fiorito 
et al. 2003), in which strong correlations of QPO frequency with 
power-law index have been observed; such states have tradi- 
tionally been noted as low-brightness, power law-dominated 

Q 2 states (see, e.g., Belloni et al. 1999), where the disk is at least 
tenuous if not undetectable altogether. 

4. Studies tracing the QPO frequency and index over time 
between state transitions show that the frequency tracks the 
index during state evolution (Homan et. al. 2001 ; Sobczak et al. 
1999, 2000) from a value of about 1.7 ± 0.1 to about 2.7 ± 0.2, 
as the QPO frequency varies from about 1 to 1 0 Hz. 

The persistence of the correlation of index with QPO fre- 
quency and its tracking with respect to time suggests that the 
underlying physical process or condition that gives rise to the 
low-frequency QPO are tied to the corona, which gives rise to 
the power-law component, and furthermore that this physical 
process or condition varies in a well-defined manner as the 
source progresses from one state to another. Moreover, the fact 
that the same characteristics of the correlation are seen in so 
many Galactic X-ray binary BH sources, which vary widely in 
both luminosity (presumably with mass accretion rate) and 


state, suggests that the physical processes controlling the in- 
dex and the low-frequency QPOs are a characteristic property 
of these sources and, by virtue of the low-high frequency 
correlations of Psaltis et al. (1999) and Belloni et al. (2002), 
may be a universal property of all accreting compact systems. 

The above data have motivated us to develop a detailed 
model of the physics of the corona surrounding a BH that 
directly predicts the behavior of the spectral index with fun- 
damental properties of the corona. The model that we have 
developed incorporates fundamental principles of fluid me- 
chanics, radiative transfer theory, and oscillatory processes. It 
identifies the origin of the QPO as a fundamental property of 
a compact coronal region near the BH and shows how the 
photon index of this corona changes as a function of mass 
accretion rate. 

It has been already shown in a number of papers (see be- 
low) that the absence of a firm surface in BHs leads to the 
development of the very strong converging flow when the 
mass accretion rate is higher than Eddington rate Me dd = 
^Edd/c 1 2 3 4 - The main observational features of the converging 
flow should be seen in high/soft phase, while the thermal 
Comptonization spectrum of the soft (presumably disk) pho- 
tons should be seen in the hard phase. 

In series of papers (Chakrabarti & Titarchuk 1995; 
Titarchuk et al. 1996, 1997; Ebisawa et al. 1996; Titarchuk & 
Zannias 1998; Shrader & Titarchuk 1998, 1999; Borozdin Q3 
et al. 1999; Laurent & Titarchuk 1999, hereafter LT99; Laurent 
& Titarchuk 2001; Titarchuk & Shrader 2002; Turolla et al. 
2002), the authors argue that the drain properties of the black Q4 
hole horizon are necessarily related to the bulk inflow in BH qs 
sources and that the spectral and timing features of this bulk 
inflow arc really detected in X-ray observations of BHs. The 
signatures of the inflow are as follows: 

1. An extended steep power law in which the photon index 
saturates to the value 2.7 ± 0.2 as the mass accretion rate 
increases; the precise value of the index is a function of the 
temperature of the flow, but its relatively high value is a result 
of inefficient photon upscattering in the converging flow due to 
photon trapping. LT99 demonstrated that the photon spectral 
index varies in a narrow range from 3 to 2.7, with mass ac- 
cretion rate (in Eddington units) increasing from 2 to 7. 
Moreover, LT99 demonstrated using Monte Carlo simulations 
that in the low/hard state when the plasma temperature is of 
order 50 keV the bulk inflow spectrum is practically identical 
to the thermal Comptonization spectrum (there is no any no- 
ticeable effect of the bulk Comptonization in the spectrum). In 
fact, the effect of the bulk Comptonization compared to the 
thermal one becomes stronger when the plasma temperature 
drops below 10 keV. The small variation of the photon index 
around 1.7 is a characteristic signature of the low/hard state that 
was pointed out in earlier work by Chakrabarti & Titarchuk 
(1995) and later confirmed by LT99. 

2. The QPO high frequency (100-300 Hz), which is in- 
versely proportional to BH mass. It is worth noting that the 
QPO frequency scaling as l/M is a generic feature of any QPO 
model in which QPO frequency is scaled with the Schwarzchild 
radius. 

Furthermore, Titarchuk et al. (1999) presented observational 
evidence that the variability of the hard X-ray radiation in NSs 
and BHs occurs in a bounded configuration. Ford & van der 
Klis (1998) found that the break frequency in the power density 
spectrum (PDS) is correlated with QPO frequency for some 
particular NS sources (4U 1728—32). Wijnands & van der Klis 
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Fig. 1 . — TL (corona) concept. This picture renders the gravitational energy release in the disk and corona and the coronal illumination by the disk soft flux in 
compact object (NS and BH). See text for details. 


(1999) later found a similar correlation in BH sources. 
Titarchuk & Osherovich (1999, hereafter T099) and Titarchuk 
et al. (1999) explained this correlation in NSs and BHs, re- 
spectively. They argued that the observed PDS is the power 
spectrum of an exponential shot, exp(-f/^), which is the 
response of the diffusion propagation to any perturbation in a 
bounded medium. The inverse of the characteristic diffusion 
time v\> = 1 jtd and the normal mode QPO frequency vq have a 
well-defined relation to the size of the bounded configuration, 
because v h oc V /L and uq oc F/Z 2 , where V is the specific 
perturbation (magneto sonic) velocity and L is a characteristic 
size of the configuration (see T099). 

Titarchuk et al. (1998, hereafter TLM98) proposed that this 
bounded configuration (cavity) surrounding compact objects 
is the transition layer (TL) that is formed as a result of dy- 
namical adjustments of a Keplerian disk to the innermost sub- 
Keplerian boundary conditions. They argued that this type 
of adjustment is a generic feature of the Keplerian flow in 
the presence of the sub-Keplerian boundary conditions near 
the central object and that it does not necessarily require the 
presence or absence of a hard surface. TLM98 concluded that 
an isothermal sub-Keplerian TL between the NS surface and 
its last Keplerian orbit forms as a result of this adjustment. The 
TL model is general and is applicable to both NS and BH 
systems. The primary problem in both NS and BH systems is 
understanding how the flow changes from pure Keplerian to 
sub-Keplerian as the radius decreases to small values. TLM98 
suggested that the discontinuities and abrupt transitions in their 
solution result from derivatives of quantities such as angular 
velocities (weak shocks). They were the first to put forth the 
possibility of TL formation to explain most observed QrOs 
in bright low-mass X-ray binaries (LMXBs). In Figure 1 we 
illustrate the main idea of the TL concept for NSs and BHs. 

In this paper we explain the general correlation between 
photon index and low frequency in the framework of the TL 
model and we give more arguments for the nature of the 
spectral phases (states) and phase transitions observed in BHs. 


The main features of the TL model are given in § 2.1. The 
formulation of the problem of low-frequency oscillations and 
the relationship with MA oscillations in the TL are also de- 
scribed in § 2.1. The coronal model and the spectral index- 
optical depth relation are given in § 2.2. The derivation of the 
low frequency-index correlation and the details of modeling 
the spectral phase transition are presented in § 2.3. We analyze 
specific QPO and spectral data in terms of the TL model in § 3. 
We discuss the signatures and the methods of the identification 
of BH and NS sources using timing and spectral character- 
istics in § 4. Our summary and conclusions follow in § 4. 

2. TRANSITION LAYER THEORY 

2.1. The Main Features of the Transition Layer Model 

TLM98 define the TL as a region confined between the 
inner sub-Keplerian disk boundary and the first Keplerian 
orbit (for the TL geometry, see Fig. 1 of this paper and Figs. 1 
of TLM98 and Titarchuk et al. 1999). The main idea at the 
basis of our investigation is that a Keplerian accretion disk is 
forced to attain sub-Keplerian rotation close to the central 
object (a neutron star or a black hole). The transition between 
the Keplerian and the sub-Keplerian flow takes place in a 
relatively narrow, shocklike region (the size of the region if of 
order the radius of the central object [NS or BH]) where 
dissipation occurs. As a consequence, the gas temperature 
increases and the disk puffs up, forming a hot corona. Then the 
corona intercepts the soft disk photons and upscatters them via 
thermal and dynamical Comptonization, finally giving rise to 
the hard power law tail. 

Tl _ . 1 'J. 1-._ _1_ J _ J _ i.1 1 

me puwci-iaw uiueX vaiuc Miuugiy ucpciiu^ uii uic euiOfim 
optical depth and temperature (see Chakrabarti & Titarchuk 
1995; Titarchuk et al. 1997; LT99). When the mass accretion 
rate in the disk increases and consequently the disk soft photon 
flux increases, the corona is drastically cooled down to a 
temperature of order 5-10 keV (in § 2.2 we present a detailed 
analysis of this effect). For such a low plasma temperature the 
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bulk motion of the converging flow is more efficient in up- 
scattering disk photons than thermal Comptonization. Fur- 
thermore, the indices saturate to the asymptotic values around 
2.7 as the mass accretion rate of the converging flow increases. 
The exact asymptotic value of the index is determined by the 
plasma temperature only (LT99). 

TLM98 evaluate the size of the TL L as a function of a 
nondimensional parameter 7 (often called the Reynolds num- 
ber), which is the inverse of the a -parameter for the accretion 
flow (Shakura & Sunyaev 1973), 

7 = M/ATTrjH — RV r /D, (1) 

where M = 4itRHm p N e V r is the accretion rate in the disk, H 
is a half-thickness of the disk, V r is a characteristic radial 
velocity at a given radius R in the disk, N e is electron (proton) 
number density, m p is the proton mass, and D is the diffusion 
coefficient. D can be defined as D = V t l t j 3 using the turbulent 
velocity and the related turbulent scale, respectively, or as 
D = D m = c 2 / a for the magnetic case, where a is the con- 
ductivity (e.g., see details of the D definition in Lang 1998). 
It is worth noting that the viscosity 77 is a function of M . 

Because 7 and the TL Thomson optical depth r 0 , which is a 
product of L = R out - R m (Rm and R out are the TL inner and 
the outer radius, respectively), electron number density N ei 
and Thomson cross section 07 (r 0 = a T N e L), are both related 
to M, one can in principle determine r 0 as a function of the 
7 -parameter. Because of the uncertainties of the disk equation 
of state and the viscosity 77 in the disk, there is still uncertainty 
in the precise determination of this Tq( 7 ) relation. We address 
this issue of the r 0 ( 7 ) determination in § 2.3. 

TLM98 determine the L( 7 ) dependence using the equation 
of motion in the disk where the radial motion is controlled by 
friction and the angular momentum exchange between adja- 
cent layers, resulting in the loss of the initial angular mo- 
mentum by accreting matter (e.g., Shakura & Sunyaev 1973). 
We define an adjustment point where a Keplerian accretion 
disk is forced to attain sub-Keplerian rotation close to the 
central object (NS or BH). In order to determine the radius of 
the adjustment point one should combine this equation of 
motion with three boundary conditions: one at the innermost 
disk radius R in and two at the adjustment radius R ou t- Namely, 
at R m the disk rotational frequency matches the rotational 
frequency of the central object (spin), and at R out there is no 
break in the rotational frequency, and its derivatives on the 
left-hand and right-hand sides are equal. Thus, for a given 
value of the 7 -parameter the rotational frequency profile u(R) 
and the outer radius of the TL are uniquely determined by the 
three boundary conditions and the equation of motion (see, for 
example, eqs. [ 8 ] and [9] of T099). 

The adjustment of the Keplerian disk to the sub-Keplerian 
inner boundary creates conditions favorable for the formation 
of a hot plasma outflow at the outer boundary of the TL, be- 
cause the Keplerian motion (if it is followed by sub-Keplerian 
motion) must pass through the super-Keplerian centrifugal 
barrier region (TLM98). 

The equation (see TLM98, eq. [11], or T099, eq. [9]) 
determines R out as a function of 7 . For given values of R out ( 7 ) 
one can find values of the Keplerian frequencies ^(Rout) (see 
eq. [1] of T099). For example, for the observed range of the 
high QPO frequencies in BH sources, from 1 to 100 Hz (see 
Fig. 2), one can find using the TLM98 solution for BH mass of 
10 M 0 that 7 varies from 2 to 20 (a variation that is related to 
the mass accretion rate M; see eq. [1]). We remind the reader 



Fig. 2. — Correlation between frequencies of QPO and noise components 
for white dwarf {diamonds , SS Cyg [Mauche 2002]; squares , VW Hyi [Woudt 
& Warner 2002]), neutron star ( open circles ), and black hole candidate {filled 
circles ) sources. Neutron star and black hole data are from Belloni et al. 
(2002). The dashed line represents the best-fit of the observed correlation (see 
Titarchuk et al. 2001b and TW02, where this correlation is predicted and 
explained using the TL model). This plot also appears in Mauche (2002). 

that in the framework of TLM98 these high QPO frequencies 
are interpreted as Keplerian frequencies i/ K at R out . 

It is worth noting that in the literature a = is treated as 
a constant and a free parameter independent of the mass ac- 
cretion rate (or the spectral state), whereas we find that a 
varies from 0.5 to 0.05 with Af, i.e., that a is a very strong 
function of M. 

Soria (1999) inferred the effective value of the viscosity a eff 
parameter as a function of the magnetic field. He showed that 
regardless of the true viscosity in the disk the a e ff parameter is 
always higher than 0.05 if the ratio of the magnetic pressure 
P m to the gas pressure P g is about 1 0% and more. Below (see 
also TW02 for details) we show, using the low frequency- 
high frequency correlation, that the inferred ratio of P m /P g = 
\/P = 0.1 is consistent with values of a =0.05-0.5 and 
Soria’s predictions of high values of a: at P m /P g £ 0.1. 

The TL model identifies the low-frequency QPO as that 
associated with the viscous MA oscillation of the bounded TL; 
this mode is common to both BHs and NSs. The correlation of 
the MA frequency with the gravitational (Keplerian) frequency 
i/K = 2nco K was derived by TO99 and Titarchuk et al. (2001b) 
for NSs, and then it was generalized by TW02 for BHs and 
white dwarfs. The MA frequency is derived as the eigen- 
frequency of the boundary-value problem resulting from an 
MHD treatment of the interaction of the disk with the magnetic 
field. The problem is solved for two limiting boundary con- 
ditions that encompass realistic possibilities. The solution 
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Fig. 3. — QPO low frequency vs. 7 -parameter (xM). 


yields a velocity identified as a mixture of the sound speed and 
the Alfven velocity. The Titarchuk et al. (2001b) treatment 
does not specify how the eigenfrequency is excited or damped. 
However, it makes clear that the QPO is a readily stimulated 
resonant frequency (see Titarchuk et al. 2002 and Titarchuk 
2003 for details of excitation of the eigenfrequencies). 

A linear relation was derived between v ma and i/ K (see Fig. 2 
and TW02), 

J'MA = CmA^K, (2) 

where Cma = 2 1/2 4tt[(1 -F /£)/(! -I- f3)]^ 2 (H/r 0Xlt ) is a pro- 
portionality coefficient that is “universal” to the extent that (3 
and HlR 0Ut remain about the same from one source to the next. 
Here (3 = P 9 /P m is the ratio of the gas pressure to the mag- 
netic pressure, the coefficient / = 1/32 and 1/327T 2 for stiff 
and free boundary conditions, respectively, and H is the half- 
width of the Keplerian part of the disk. 

The identification of i^igh and v ]ow with i/ K and vma> re- 
spectively, leads TW02 to the determination of H/R 0 ut = 
1.5 xl0~ 2 and f3 = 10 using the observed value of Cma = 
0.08. On the other hand, uk = ^high at the adjustment radius 
R oni as a function of the 7-parameter can be readily obtained 
using equation (9) of T099 [see above for details of ^(7) 
determination]. Thus, as a result of the identification of v j ow 
with ^ma> we obtain the relation z/j ow = 174 A = Cma^ as a 
function of 7, namely, 

i/ low = 0.08 i/ k (7)- (3) 


o 



equations (1) and (9) of T099, keeping in mind that i/ high = 
z/ K . Because the radius of interest, which in our case is the 
outer TL boundary radius i? out , is related to i/ K > one can relate 
i? ou t to I'low using formula (3). We present this function in 
Figure 4. 

The spectral index-frequency correlation can be derived if 
we find a spectral index-TL optical depth correlation. Because 
7 is a function of the mass accretion rate M and the optical 
depth tq is a function of M , we are led to a correlation of y low 
and ro and ultimately to the index-frequency correlation. It is 
worth noting that the 7 and ro dependences on M imply that a 
correlation between these quantities exists but do not imply 
that there is a linear correlation between 7 and To, because the 
7 - parameter also depends on the accretion flow viscosity r), 
which is also a function of M. 


In Figure 3 we present the inferred z^ ow as a function of 7. 
We derive this relation using equation (9) in T099 with 
the assumptions that a K = m(y 0 /363 Hz) = 1 (see eq. [7] of 
T099 for definition ajQ and m = M/M , 0 = 10. 

It is evident that the inferred z^high and *7 0W are inversely 
proportional to m because 

tfrugh = v*. °c (4) 

i/low (5) 

where L is the TL characteristic size, which is proportional 
to m. 

We calculate the outer radius of the TL, r out = R ou t/Rs, 
as a function of i/j ow using equation (2) or (3) above and 


2.2. Spectral Index-Optical Depth Relation: Corona Model 

TLM98 offered the TL model adjustment of the Keplerian 
accretion flow (disk) to the sub-Keplerian inner boundaiy 
condition, which provides a compact corona for Comptoni- 
zation upscattering of the illuminating accretion disk photons. 
The essentially sub-Keplerian rotation of NSs in LMXBs is a 
well-established observational fact due to the direct mea- 

aUiLiuviiid ui apiua a icvicw uy van uci ivui -4UUU JLUI 

details). Direct BH spin measurements are still unavailable. 
Because the bulk motion spectral and timing model for a 
Schwarzschild black hole is probably consistent with obser- 
vations (see Laurent & Titarchuk 2001 for details), one can 
conclude that a fast-rotating Kerr BH is ruled out by the 
observations. In fact, there is no space between the inner edge 
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of the disk and the horizon in the fast-rotating Kerr BH to 
upscatter the disk soft photons by the bulk motion to the en- 
ergies of order 500 keV that are detected in high-energy 
observations (Grove et al. 1998). Thus we can assume that 
BHs may also undergo sub-Keplerian rotation. 

It is likely that this adjustment to sub-Keplerian flow is not 
smooth and that near the adjustment radius the strong or weak 
shocks can be formed. The adjustment radius (shock) region 
can be treated as a potential wall at which the accreting matter 
releases its gravitational energy in a geometrically thin target. 
As a consequence of this the shock plasma temperature is 
much higher than that in the surrounding regions of the 
Keplerian disk. The shock region gets puffed up oscillating 
with frequencies close to (Titarchuk 2003). Additional 
oscillations can appear (at least in the NS case) because of the 
rotational configuration above the disk, which is a magneto- 
sphere in the NS case (see details in Titarchuk et al. 1999). 
The shock formation leads to the formation of the hot inner 
bounded region. The hard photons in the shocked region il- 
luminate the surrounding material, evaporating some part of 
the disk. It is obvious that the evaporated fraction depends 
on the mass accretion rate in the disk (see Chakrabarti & 
Titarchuk 1995). Some small fraction of the hard photons (at 
a maximum of 6%) can be reflected by the cold parts of the 
disk (Chakrabarti & Titarchuk 1995; see also Basko et al. 
1974). 

Let us assume the column density of this region y 0 is a few 
grams or that the Thomson optical thickness tq is a few (see 
Shakura & Sunyaev 1973 for the thickness estimate of the 
innermost part of the disk). The rate of energy release at this 
region g cor is a few percent of the Eddington luminosity since 
the adjustment radius is located within 4-100 Schwarzschild 
radii (see Fig. 4). The heating of gas due to the gravitational 
energy release should be balanced by the photon emission. For 
high plasma temperature, Comptonization and free-free emis- 
sion are the main cooling channels, and the heating of electrons 
is presumably due to the Coulomb collisions with protons (see 
Fig. 1 for a picture of the disk corona). 

Under such physical conditions the energy balance can be 
written as (see also Zel’dovich & Shakura 1969; TLM98) 

Qcor/^O — C eomp e(r)T e /f(T e ) + C e T l J 2 p, (6) 

where r is the optical depth within the shock, e(r) is a dis- 
tribution function for the radiative energy distribution, f(T e ) = 
1 4- 2.5 (kT e /mc 2 ) is the relativistic correction factor, T e is the 
electron temperature in kelvins, and C c0 mp = 20.2 cm s' 1 K" 1 
and Cff = 2.6 x 10 20 ergs K _1/2 cm s” 1 are dimensional con- 
stants. In this formula we neglect gas heating due to recoil 
effect. 

The distribution e(r) is obtained from the solution of the 
diffusion equation (Zel’dovich & Shakura 1969), 

dh/dr 2 = -3(0 tot /c)/ro, (7) 

subject to two boundary conditions. Here we assume that 
the region of the gravitational energy release (corona) is a 
spherical shell surrounding the central object and the total flux 
in the corona is a sum Q to t = Qcox 4- Qd of the gravitational 
energy release Q cor and the illumination flux from outside of 
the corona (disk) Q d . The inner and outer boundary condition 
for a BH corona is that there is no scattered radiation from 
outside of the corona. Then in the Eddington approximation 


(see Sobolev 1975) the boundary conditions can be written as 
follows: 

(de/dr - 3 € /2)| t=0 = {de/dr + 3e/2)| T=To = 0, (8) 

where r — 0 and r = to are at the inner and outer boundaries, 
respectively. 

For an NS corona one should assume that inside of the inner 
corona boundary the energy density e = constant, namely, that 
the radiation flux emitted toward the central object at some 
point of the inner boundary returns back (reflected by the 
NS surface); this condition is equivalent to the reflection 
condition: 

de/dr\ T=0 = 0. (9) 

Below we present a general formula of the coronal tempera- 
ture that combines the BH and NS cases. 

The solution of equations (7) and (8) provides us with the 
distribution function for the energy density in the BH corona: 

t(r) = (2,o./c){ 1 + (3/2 )t 0 [t/t 0 - (t/t 0 ) 2 ]}. (10) 

We neglect the dependence of Q tot on r to derive this formula. 
Furthermore, we take a representative value of (e) for the 
maximum value in the corona that is very close to the mean 
value of 6 mean , i.e., 

(c) = 6 max = (e to t/c)[l+(3/8)r 0 ], (11) 

in order to approximate the coronal temperature. In Figure 5 
we present the results of calculations of the temperatures kT e 
(keV) and the energy spectral indices of the Comptonization 
spectrum a (photon index T = a 4- 1) as a function of the 
optical depth of the shell r 0 . The different curves of r e (r 0 ) are 
related to the different ratios of QJQc or- The calculations have 
been made for values of p — 10 -6 g cm" 3 and Q cor = 10 22 ergs 
cm -2 s~\ which are characteristic values of the density and lu- 
minosity for the standard disk model (Shakura & Sunyaev 1 973). 
In fact, the calculation results weakly depend on p and Q C0T if 

(p/10~ 6 g cm“ 3 )/(£> cor /10 22 ergs cm -2 s _l ) S 1 

is of order a few or less. In this case the equation for the 
temperature is simplified 

6.73 x 10- 2 t 0 (1 + 3t 0 /8)(7;/10 8 K)(l + Q d IQ, 0 ,)lf{T e ) = 1, 

( 12 ) 

which has the solution 

Tel 10 8 K= 1.5[t 0 (1 + 3ro/8)(l + &/&or) - 0.62]" 1 . (13) 

It is easy to show that the case for which Qd/Qco r = 1 is 
identical to the NS case (the inner reflective boundary) when 
the disk illumination is neglected. One can see from equa- 
tion (13) that T e is a function of To if Qd/Qcor 1. It is also 
worth noting that T e ro(\ 4- 2ro/S)/f(T e ) is insensitive to the 
total luminosity 0tot if Qd/Qcor <£ 1 (eq. [12]). 

The spectral index a as a function of the product of T e r 0 is 
insensitive to £>tot too for r 0 ^ 1 because the Comptonization 
parameter 

y = 4 kT e /m e c 2 max (r 0 , Tq) 
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Fig. 5. — Energy spectral index vs. the TL optical depth (thermal Comp- 
tonization case). 

and 


a = -3/2 + \/(9/4 + TJy) 

(see Rybicki & Lightman 1979; Sunyaev & Titarchuk 1980). 
This inferred property of the Comptonization spectra repro- 
duces the observed independence of a on the bolometric lu- 
qs minosity (see, e.g., Tanaka & Lewin 1995) when a varies 
within the range 0.6 ± 0.2 (T = 1.6 ± 0.2). This defines the 
so-called low/hard state in BHs. 

Below we show that the observed plateau of the index- 
frequency correlation at low values of Vi < 1 Hz is the result of 
this behavior because the lower frequency values are related to 
the low mass accretion rates m — M/M^dd < 1 and r 0 is of 
order 1 when m < 1. 

We calculate the exact values of the spectral indices a for a 
given T e and r 0 using the relativistic formulae developed by 
Titarchuk & Lybarskij (1995). See their formulae (17) (for 
plane geometry; note that their r 0 is our r 0 /2) and (24). In 
Figure 5 we present the spectral indices as a function of r 0 for 
BH and NS cases when the disk illumination is negligibly 
small with respect to the coronal energy release (a condition 
that can be a model for the iow/hard state in these systems). 
The plasma temperature values of 20-150 keV for r 0 = 1-3 
and a = 0.6 ± 0.2 (for ratio = Qd/Qcor = 0) are typical val- 
ues of these quantities for the low/hard state in BHs. 

When Q d is comparable to g cor , the cooling becomes more 
efficient because of Comptonization and free-free processes 



Fig. 6. — Photon spectral index vs. the TL optical depth (general case). 


and therefore T e unavoidably decreases (see T e values for 
Qd/Qcor = 1-4 in Fig. 5). 

This illumination effect and mass accretion rate increase 
may explain the hard-soft transition when the Compton tem- 
perature drops substantially with increase of the soft (disk) 
photon flux (for illustration, see Fig. 1). In Figure 5 we show 
that the temperature drops from 50-60 keV for ratio = 0 to 7- 
1 0 keV for ratio = 4. The temperature may also drop because 
of an increase in the optical depth (presumably because m 
increases) even when ratio = 0. 

2.3. Modeling of Spectral Phase Transition from the Hard 
State to the Soft State 

LT99 studied the Comptonization of the soft radiation in the 
converging inflow (Cl) into a BH using Monte Carlo simu- 
lations. The full relativistic treatment has been implemented to 
reproduce the spectra. LT99 show that spectrum of the soft 
state of BHs can be described as the sum of a thermal (disk) 
component and the convolution of some fraction of this 
component with the Cl upscattering spread (Green’s) function. 
The latter boosted photon component is seen as an extended 
power law at energies much higher than the characteristic 
energy of the soft photons. LT99 demonstrate the stability of 
the power-law index (the photon index T = 2.8 ± 0.1) over a 
wide range of the plasma temperature 0-10 keV and mass 
accretion rates (higher than 2 in Eddington units) due to 
upscattering and photon trapping in the CI. The spectrum is 
practically the same as that produced by standard thermal 
Comptonization when the CI plasma temperature is of order 
50 keV (the typical ones for the BH hard state) and the photon 
index T is around 1.7. 

LT99 also demonstrate that the change of the spectral 
shapes from the soft state to the hard state is clearly related to 
the temperature and optical depth of the bulk inflow. We 
combine our results of calculations of index and Compton 
cloud temperatures for thermal Comptonization case (see § 2.2 
and Fig. 5) and LT99’s results of the spectral calculations (see 
Table 2 of LT99 for the spectral index values) to evaluate the 
power-law index as a function of the optical depth of the 
Compton cloud r 0 . We use the values of the index and of 
the temperature for Q = 0 (see Fig. 5) when r 0 varies from 1 
to 3 to describe the index behavior and the plasma temperature 
variation in the hard state. Figure 6 presents the photon index T 
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as a function of tq. This function F(r 0 ) can be fitted analyti- 
cally as follows: 

_ f C° + CiTo + c 2 t 2 0 + c 3 rl 1< r 0 <5.17, 

(T0) “ l b 0 - VVn To > 5 . 17 , (14) 

where Co — 0.61, cj — 0.57, C 2 = 1.485 x 10~ 2 , C 3 = —8.5 x 
10- 3 , b 0 = 2.81, b 1 = 1.0 xlO 3 , and r, =0.5. 

In order to derive the index dependence T on one 
should relate the parameter 7 to r 0 and then combine the 
two derived relationships T versus To and ^i ow versus 7 (see 
Figs. 3-6). The optical depth r 0 and 7 -parameter are functions 
of M. To fit the data we assume that r 0 — Aj 6 . We have three 
Q9 free parameters, BH mass m — M/M { ©, 6, and A , to fit the 
model to the data if we assume that the Cl is essentially cold 
(the Cl temperature kT c \ ^ 5 keV). This relationship is all that 
is needed to fit the observed T-i/i ow correlations. 

3. DATA AND ANALYSIS 

We analyze data points for photon index versus low fre- 
quency correlation obtained from Trudolyubov et al. (1999) 
and V03 for GRS 1915+105 and Sobczak et al. (1999, 2000) 
for XTE J1 550-564. Also, we make use of data for GRS 
1915+105 by Fiorito et al. (2003), which is a detailed study 
of QPO and spectral features of the steady low/hard states 
of GRS 1915+105 for representative 1996 and 1997 obser- 
vations. Fiorito et al. computed the power spectra of combined 
binned and event mode data for PCA channels 0-255 and 
sampled at 0.25-64 Hz for 4 s intervals. The spectra were 
averaged over the entire observation duration. The PDS spectra 
were fitted with a combination of Lorentzian lines and ex- 
hibited a break frequency, a central QPO frequency, and in 
some case the presence of the first harmonic. The background 
subtracted energy of each observation ID was obtained and 
time-averaged over the entire observation duration using stan- 
dard 2 PCA and HEXTE data where available (both clusters). 
The XSPEC model used to fit the energy spectra consisted 
of a multitemperature thermal disk component (DISKBB), a 
simple power law, and a small iron line component. For the 
1996 data (ObsSID 10408), it was found that the disk com- 
ponent could be neglected without affecting the fit. Accept- 
able x 2 values were obtained in all cases. 

3.1. Correlations between QPO Frequency and 
Spectral Parameters 

The first part of the Fiorito et al. (2003) study was to look 
for correlations of different spectral components, presumably 
thermal disk components, with QPO frequency to try to de- 
termine the reason for the nonunique total flux correlations 
previously observed. In the case of ObsID 20402 (1997), we 
did observe the expected positive correlation of QPO fre- 
quency with disk flux, in agreement with previous analyses of 
qio Muno et al. (1999) and Markwardt et al. (1999), which were 
done in other states. However, we observed an anticorrelation 
of frequency with disk flux for ObsID 10408 (1996) and no 
correlation of frequency with disk flux for ObsID 10258 
(1996). Furthermore, no correlation of QPO frequency with 
disk temperature is seen for any of the ObsIDS we studied. 
These results are at variance with the often-quoted result that 
disk parameters are closely associated with QPO frequency. 
Remillard et al. (2002b) show a comparison of low QPO 
frequency and the apparent disk flux for the source XTE 
J1 550—564 during its 1998-1999 outburst. The QPOs are 



0.2 0.5 1 2 5 

QPO low frequency (Hz) 


Fig. 7. — Power-law index vs. QPO centroid frequency for the plateau 
observations of GRS 1915+105 from V03, along with a fit using the TL model 
with m = 12 and to — 7 1 25 . 

reported with differences in phase lags. Again, it is clear that 
the frequency-disk flux correlation observed for this source is 
only observed at lower (u < 5 Hz) QPO frequencies and 
fluxes than usually observed for GRS 1915+105. 

However, robust correlation between index and the central 
frequency of low-frequency QPOs is also seen in a wide variety 
of BH sources radiating in different states and during transition 
between states. The QPO frequencies in such studies are ob- 
served in the PDS obtained by integrating over time intervals 
ranging from the complete time of the observation to segmented 
time bins where the luminosity remains constant (see, e.g., 
Markwardt et al. 1999). Strong frequency-index correlations 
are observed by Kalemci ( 2002 ) in outburst decay (presumably 
state transitions) of a number of sources and by V03 in a 
comprehensive study that focuses its attention on the behavior 
of the power-law index and QPO frequency in several BH 
candidate sources. We will use the data of V03, our studies, and 
those of other studies as well in our analysis (see § 3.2 below). 

Observation of low-term tracking of the QPO frequency 
with the variations of the photon index has been observed in at 
least two BH sources. Tracking is seen in the data of Rossi 
(2003) and Homan (2001) in observations of XTE J1 560-500 
taken over a 30 day period. In this period the source appar- 
ently transits from the low/hard state (T ~ 1 .7) to the inter- 
mediate or soft state, where the index apparently remains close 
to the constant at a value (T ~ 2.7). In this transition the QPO 
frequency varies from about 1 to 10 Hz. Similar tracking is 
observed in XTE J1 550-564 by Sobzcak et al. (1999, 2000). 
However, the index and QPO frequency does not track with 
the total X-ray flux (V03, Fig. 7). 

3.2. Comparison of the Inferred Index-Frequency Correlation 
with the Data: Low/Hard-High/Soft State Transitions 

3.2.1. GRS 1915+105 

In Figure 7 we present the results of fit to the data of V03 
for the plateau observations of the microquasar GRS 
1915+105 using the TL model. The best-fit parameters are 
m = 12, 6 = 1.25, and A = 1 . We note that the value m = 12 
is consistent with the BH mass evaluation m = 13.3±4 
obtained by Greiner et al. ( 2001 ) using the IR spectra and by 
Shrader & Titarchuk (2003; see also Borozdin et al. 1999) 
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Fig. 8. — Power-law index vs. QPO centroid frequency for the observations 
of classes /3 and u {black points) and a and u {red points ; observations 15 and 
16) of GRS 1915+105 from V03, along with a fit using the TL model with 
m = 12 and r 0 = 7 125 . Values for the plateau observations (see previous 
figure) are plotted for comparison {blue points). 


using the X-ray spectra. The observable index saturations at 
low and high values of v L are nicely reproduced by the model. 
The low- and high-frequency plateau regions of the data and 
fit are signatures of the two spectral phases that are explained 
by two different regimes of Comptonization (upscattering): 
( 1 ) bulk flow Comptonization in the soft state (saturation at 
T ~ 2.8) and (2) thermal Comptonization in the hard state (the 
photon index tends to level at F ~ 1.7). 

A comparison of the index-frequency correlation of plateau 
observations with V03’s (3 and v data (Belloni 1999 classifi- 
es cation scheme) and the model fits are plotted in Figure 8 . The 
only difference in the data is different saturation values at high 
v L , which can be readily explained by a change of plasma 
temperature in the converging inflow kT ci . For the upper value 
of T = 2.81 the temperature kT ci = 5 keV, and for F — 2.7 the 
temperature kT c \ = 9 keV. The best-fit parameters for the latter 
data are the same as for the plateau observations, except that 
Q 12 in the model function T (r 0 ) (see eq. [14]) the coefficients b o, 
b\ 9 and r* are replaced by bo = 2.7, b\ = 22.4, and r* = 0.8 
for r 0 > 2.47. 

In Figure 9 we present a plot of power-law index versus 
QPO centroid frequency for the observations of class a and v 
of GRS 1915+105 from V03 {black points; observation 18). 
Blue points correspond to the values for observations by 
Fiorito et al. (2003) for ObsIDs 20402 and 10258. Magenta 
points correspond to positions where the ObsID 20402 data 
points are shifted to one-half the measured QPO low-frequency 
values. This suggests that the latter data points are related to the 
second harmonics of the v L frequencies. It is well known that 
data for NS QPOs are a mixture of the first and the second 
harmonics of HBO frequencies (see, e.g., van der Klis 2000). 
Red points correspond to positions where there should be 
points with half of the frequencies for V03’s observation 1 8. A 
theoretical curve {blue solid line) obtained using F(r 0 ) (see 
eq. [14]; for coefficients bo — 2.7 , b\ = 22.4, and r* = 0.8 for 
ro > 2.47) form = 12 and To = 7 1 5 fits Fiorito’s ObsID 10258 
data points. 

It is worth noting that relation between 7 and r 0 is slightly 
different from what we obtained for the best fits presented in 
Figure 8 . Namely, the index 7 increases from 1.25 to 1.5. This 



QPO low frequency (Hz) 


Fig. 9. — Power-law index vs. QPO centroid frequency for die observations 
of classes a and v of GRS 1915+105 from V03 {black points', observation 18). 
Blue points correspond to the values for observations by Fiorito et al (2003; 
see text). Magenta points correspond to positions where there should be points 
with half of those frequencies. Red points correspond to positions where there 
should be points with half of the frequencies for observation 18. The blue 
solid curve is for m — 12 and r 0 = 7 1 5 . 


means that for the same value of 7 the optical depth to is 
slightly higher than that in the state presented in Figure 8 . This 
change can be explained by a slightly higher accumulation of 
plasma in the TL. For example, if we fix the disk viscosity 7 7 , 
then the 7 -parameter is determined by M only. But the mass 
accretion rate M depends on the product of number density N e 
and the flow radial velocity V r , i.e., M oc N e V r . Thus, for a 
given value of NV r the optical depth to oc N e varies as V r . 
Smaller values of V r corresponds to the higher values of N e and 
consequently to higher values of t 0 ; namely, higher plasma 
accumulation in the TL can be related to smaller velocities V r . 

It is interesting that the suggested half-frequency points 
{magenta and red points) are located in a narrow corridor 
between these inferred curves (presented in Figs. 8 and 9, 
respectively). 

3.2.2. XTE J1 550-564 

Figure 10 shows a plot of power-law index versus QPO 
frequency for XTE J1 550—564. All data points presented in 
Figures 10 and 11 (see also Figs. 6 and 8 of V03) were 
obtained by Sobczak et al. (1999, 2000) and Remillard et al. 
( 2002 a, 2002 b). 

The data again show a flattening out or saturation of the 
index at a value F ~ 2.8. It is clearly seen from Figure 8 of V03 
that most of the points presented there can be fitted by a smooth 
line except for a few points on the right that we temporarily 
exclude for our consideration. The shift transformation v’ L = 
Ylvij 10 of the GRS 1915+105 curve F(i/) (see Fig. 7) into 
T(i/) produces a fit for the XTE 1550-564 data. This means 
that if two sources are in the same accretion regime [i.e., when 
r 0 ( 7 ) is the same], their relative BH masses can measured 
by ilie simple shift transformation v : L — m\u L /m 2 . In our case 
we use m \ = 12 for GRS 1915+105 and m \ = 10 for XTE Q13 
1550-564. These values are very close to the values that have 
been obtained by Shrader & Titarchuk (2003) for these sources 
using X-ray spectroscopic methods. Thus the TL model pro- 
vides an independent way to estimate the mass of one or more 
BHs in the same accretion state. 
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Fig. 10. — Comparison of the observed ( points ) and the theoretical corre- 
lations ( curves ) of photon index vs. QPO low frequency between GRS 
1915+105 (observations by V03) and XTE J1 550-564 (observations by 
Sobczak et al. 1999, 2000; Remillard et al. 2002a, 2002b; see also Figs. 6 and 
8 in V03). Black points and curve are for XTE J1 550-564; red points and 
curve are for GRS 1915+105. The XTE J1 550- 564 curve is produced by 
sliding the GRS 1915+105 curve along the frequency axis with factor 12/10 
(see text for details). 

A comparison of all data points of correlating photon index 
versus QPO low frequency for XTE 1550-564 ( black points) 
with the inferred correlation is shown in Figure 1 1 . The sat- 
uration index value in the theoretical curve is related to the 
specific value of the plasma temperature of the converging 
inflow kT cl) which changes from 5 to 20 keV from the top 
curve to the bottom curve, respectively. Slightly different 
cooling regimes of the site of the converging bulk inflow can 
explain the different saturation levels of the index observed in 
XTE 1550-564. The temperature of the converging inflow is 
determined by the disk illumination geometry and by the in- 
trinsic heating of the flow by pairs produced very close to the 
BH event horizon (see Laurent & Titarchuk 2004). All these 



QPO low frwquoncy (Hz) 

Fig. 1 1 . — Comparison of the observed correlation of photon index vs. QPO 
low frequency for XTE 1550-564 by Sobczak et al. (1999, 2000) and 
Remillard et al. (2002a, 2002b) (black points) with the inferred correlation. 
The saturation index value in the theoretical curve is related to specific value 
of the plasma temperature of the converging inflow kT ci> which changes from 5 
to 20 keV from the top curve to bottom. 


Fig. 12. — Comparison of the observed correlation of photon index vs. QPO 
low frequency for 4U 1630-47 (Trudolyubov et al. 1999; Tomsick & Kaaret 
2000; Kalemci 2002) with the inferred correlation. The saturation index value 
in the theoretical curve is related to specific value of the plasma temperature of 
the converging inflow kT C i, which changes from 12 to 20 keV from top to 
bottom. 


effects may be manifested by the observed index-frequency 
correlation at high values of v L . 

3.2.3. 4U 1630-47 

Figure 12 shows a plot of index versus QPO frequency for 
4U 1630-47 (Trudolyubov et al. 1999; Tomsick & Kaaret 
2000; Kalemci 2002). The data again show a similar behavior 
to that of XTE J1550 — 564 with an index saturation at T = 
2.45 and F = 2.65. One can see two phases with the index 
level at T — 1.7 ± 0.1 and T — 2.6 ± 0.1 and the transition 
between them that is much sharper than that observed in the 
other sources discussed above. The best-fit curve T(u L ) ( solid 
line) is obtained using to = 1 .557° 3 f° r 7 < 7 * — 3.6, to = 
1 . 557 T°* 55 7 a85 for 7 > 7 *, and the F(t 0 ) function (see for- 
mula [14]) for the best-fit BH mass m = 16, which is close 
to the estimate made by Borozdin et al. (1999) using the X-ray 
spectroscopic method. The solid and dashed curves are for 
saturation index values T = 2.65 (related to kT c i — 12 keV) 
and T = 2.45 (related to kT C] = 20 keV), respectively. Co- 
efficients 6 0 — 2.45 and 2.65, b\ — 53.7 and 44.3, and t + = 
0.55 and 0.65 for t 0 > 2.2 and t 0 > 3.42, for the solid and 
dashed curves, respectively. 

Our explanation for the observed anomalously sharp tran- 
sition is as follows. Because the sharp transition QPO fre- 
quencies do not change we can speculate that the sharp TL 
size also does not change. But if the TL size is almost the 
same, then the 7 -parameter is almost the same because the TL 
size depends on the 7 -parameter only (see § 2.1). But the 
spectral index is a function of t 0 and the TL plasma temper- 
ature. If the index rises drastically this means that upscattering 
is also drastically suppressed during the transition. Then we 
have only one possibility, i.e., that we are observing a sharp 
switching from the thermal Comptonization to the bulk mo- 
tion Comptonization regimes. In the end of this transition we 
really see a signature of the bulk inflow in terms of the index 
saturation. The bulk motion Comptonization is dominant 
when the plasma temperature drops significantly to the value 
of order 10 keV. This can happen when a strong cooling 
emerges either because the optical depth tq increases because 








No. 2, 2004 


SPECTRAL INDEX-FREQUENCY CORRELATION 


7/22/04 (V61 2/599 14) 11 


of the accumulation of matter in the TL or because an addi- 
tional source of cooling such as a soft photon supply from the 
disk appears (see investigation of these cases in § 2.2). The 
latter case can be realized if the corona intercepts more disk 
soft photons because of a puffing up of the disk or a build up 
of disk beneath the corona. 

4. DISCUSSION AND SUMMARY OF THE RESULTS 

We have developed a model of accretion onto a BH that 
greatly simplifies and reclassifies the plethora of observational 
“states” assigned to categorize the X-ray observations of 
variable BH candidates into two generic phases (states): 

1. A soft state, in which we see the effect of the BH as a 
“drain.” Bulk inflow upscattering and its quenching effect on 
the Comptonization of disk photons dominates the behavior of 
the BH spectrum. The power-law spectrum is steep in this sit- 
uation. The observed high-energy photons are emitted from 
a compact region, where soft energy photons of the disk are 
upscattered by bulk matter inflow, forming the steep power law 
with photon index around 2.8, and low QPO frequencies (above 
10 Hz) and high QPO frequencies (of the order of 100 Hz) are 
observed. In terms of the relativistic particle acceleration in the 
region of the corona and its effect on upscattering the soft disk 
photons, the soft-state spectrum is a result of the first-order 
Fermi acceleration with respect to Vic. 

We would like to emphasize that bulk inflow is present in a 
BH when the mass accretion is high, but not in an NS , where 
the presence of the firm surface leads to a high radiation 
pressure, which eventually stops the accretion. The bulk inflow 
and all its spectral features are absent in NSs; in particular, the 
saturation of the index T ~ 2.8 with respect to QPO frequency, 
which is directly related to the optical depth and mass accretion 
rate, is observed in the soft spectral state of only BHs and is 
therefore a particular signature of a BH. 

2. A hard state, which is comparatively starved for accre- 
tion. The hard phase (state) is related to an extended thermal 
Compton scattering cloud (cavity) characterized by a photon 
index F ~ 1 .7 and the presence of low QPO frequencies (below 
1 Hz). The low/hard spectrum is the result of a Fermi accel- 
eration of the second order with respect to {Vic) 1 2 3 4 5 . The effect of 
the first order on Vic is smeared out by the quasi symmetry of 
the particular dynamic, predominantly thermal motion of the 
Compton cloud plasma. 

The transition layer Comptonization model does the 
following: 

1 . It specifies the precise behavior of the power-law index 
for each state. 

2. It identifies the origin and nature of low-frequency and 
high-frequency QPOs. 

3. It predicts and explains the observed relationship be- 
tween the QPOs and their dependence on mass accretion rate 
and spectral photon index for each of the phases. 

4. It predicts and explains the observed robust correlation of 
low-frequency QPOs with the spectral index observed in BH 
candidates. 

5. It predicts and explains the observed correlation between 
low- and high-frequency QPOs in BHs, NSs, and white dwarfs. 
We interpret the correlation between low frequency V\ ow and 
power-law photon index T in QPOs extensively investigated 
by Sobczak et al. (1999, 2000), Remillard (2000b), Kalemci 
(2002), and V03 for a variety of BH sources and states, and 
also by Homan et al. (2001) and Fiorito et al. (2003). The 


observed correlation strongly constrains theoretical models 
and provides clues to understanding the nature of the QPO 
phenomena, particularly in BHs and in compact objects in 
general. 

The key assumptions and limitations of the model are as 
follows: 

1. A sub-Keplerian rotation of central object (NS and BH), 
which is an observationally well established fact in NSs of 
LMXBs but still needs verification for BHs. However, the 
existence of an adjustment region, i.e., the TL, is a quite rea- 
sonable assumption for all but very rapidly rotating BHs. Un- 
fortunately, up to now there has been no direct measurement of 
a BH spin. 

2. Our fit to the data uses the model dependence of the 
index on the TL optical depth t 0 and the model dependence of 
the QPO frequencies on the 7 -parameter. It is clear from the 
model that r 0 and 7 are correlated, but within the TLM98 
model we cannot determine the exact 7 - 0 ( 7 ) dependence be- 
cause of the uncertainty of the disk equation of the state and 
the viscosity in the disk. Therefore we have determined the 
ro( 7 ) relation by fitting the observed data, i.e., the correlation 
between the QPO low frequency and the photon spectral in- 
dex. As a result of our analysis of the observed index-QPO 
correlation we conclude that the function r 0 (j) is almost the 
same for GRS 1915+105 and XTE J1550-564 but that it is 
drastically different for 4U 1630-47. For this source the in- 
ferred tq shows a very strong dependence of 7 on the QPO 
frequencies produced during a sharp spectral state transition. 
We speculate that this peculiar behavior of the index versus 
QPO frequency and consequently the inferred r 0 versus 7 is a 
result of the plasma accumulation in the TL. It is possible that 
a similar index-frequency correlation can be detected in other 
sources. 

3. The model fits to the data give us the relatively low 
values of 7 = 2-20 (corresponding to viscosity parameter 
a = 7 _1 from to 0.05 to 0.5). Such low values of 7 can be 
understood in the framework of MHD treatment of the disk 
viscosity. In fact, Soria (1999) inferred a of order 0.1 (re- 
gardless of the true disk viscosity) and higher if in the disk a 
magnetic pressure is 10 % of the gas pressure or higher. 

4. Within the TLM98 model one can estimate the absolute 
normalization of the magnetoacoustic (QPO low) frequencies 
(see, e.g., Titarchuk et al. 2001b), but a more precise normal- 
ization can only be obtained using the observed low-high fre- 
quency correlation (see Fig. 2). 

5. In some very soft (or extended power law) regimes the 
positive correlation of frequency and index appears to be un- 
bounded, i.e., there is no indication of saturation of the index 
value (Fiorito et al. 2003). Explanation for this effect is out of 
the scope of the present TL model. We can only suggest here 
that the extended power law regime occurs when strong out- 
flows obscure the bulk inflow and the relatively cold outflow 
from winds downscatters the emerging high-energy photons 
and softens the observed spectrum (we will provide the details 
of this picture elsewhere). 

6 . V03’s Figure 6 plots the index-frequency correlation for 
the two sources XTE J1550— 564 and GRO J1 655 — 40 to- 
gether. The behavior of the later sources is quite distinct from 
all others analyzed: it shows a reverse or negative correlation 
between photon index and frequency. This is the only case 
known to exhibit this type of behavior. We propose that at a 
high accretion rate regime the pair-production heating adds to 
and can dominate the heating in the cavity, causing a transition 
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back to the thermal Comptonization phase. In this phase a 
reverse correlation between low QPO frequency can result and 
is our explanation of the single case where this phenomena is 
observed. Further theoretical investigation of this case will be 
addressed in future studies. 

The main results of this paper are as follows: 

1 . We find that the observed low frequency-index correla- 
tion is a natural consequence of an adjustment of the Keplerian 
disk flow to the innermost sub-Keplerian boundary conditions 
near the central object. This ultimately leads to the formation of 
the sub-Keplerian transition layer between the adjustment ra- 
dius and the innermost boundary (the horizon for BH). 

2. In the framework of the TL model, i^igh is related to the 
gravitational frequency at the outer (adjustment) radius v g « 
and i/jcw is related to the magnetoacoustic oscillation frequency 
i'ma- Using a relation between ^ma and the mass accretion rate 
and the photon index T and the mass accretion rate, we infer a 
correlation between z/ MA and the spectral index T. 

3. The identification of ^ low with allows us to make a 
comparison of the theoretically predicted correlation with the 
observed correlation. For this identification we use the single- 
temperature plasma assumption. We apply the plasma tem- 
perature obtained from the Comptonization spectra (electron 
temperatures) for calculations of magnetoacoustic frequen- 
cies, which strongly depends on the proton temperature. The 
single-temperature assumption is really consistent with the 
data. 

4. We present strong arguments that in BHs two particular 
distinct phases occur, where one of them, the steep power law 
phase, is a real BH signature. 

5. We found that the hard phase (state) is related to an 
extended Compton cloud (cavity) characterized by a photon 


index of around 1.7 and low QPO frequencies below 1 Hz. 
This is the regime where thermal Comptonization dominates 
the upscattering of soft disk photons and the spectral shape 
(index) is almost independent of mass accretion rate. 

6. We find that the soft phase (state) is related to the very 
compact region where soft energy photons of the disk are 
upscattered, forming a steep power law spectrum with the 
photon index saturating around 2.8. This is the regime where 
bulk flow Comptonization dominates and the effect of an in- 
crease in the mass accretion is offset by the effect of photon 
trapping in the converging flow into the BH. 

7. In the context of the distinction between NS and BH 
sources, it is worth pointing to a certain situation that can occur 
in low/hard states. The source should be an NS if the soft 
blackbody-like component with a color temperature of order 
keV is present in the spectrum (Torrejon et al. 2004). In con- 
trast, such a high temperature is necessarily related to the high 
disk luminosity in the BH case. 

8. We offer a new method of BH mass estimate using the 
index-frequency correlations. Namely, if the theoretical curve 
of the index-frequency dependence T(t/ L ) related to the pa- 
rameter m \ fits the data for a given source, then the simple slide 
of the frequency axis v' L = {rn\jmi)v L with respect to v L may 
allow us to obtain the mass m 2 by the fit of T(z/[) to the ob- 
servable correlation for another source. 


We acknowledge Tomaso Belloni for kindly supplying us 
with the data published in V03. We gratefully appreciate the 
referee’s suggestions for improving and clarifying our results 
and the limitations of our model. 
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